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Abstract —The millimeter wave (mmWave) band, a prime 
candidate for 5G cellular networks, seems attractive for wireless 
energy harvesting since it will feature large antenna arrays and 
extremely dense base station (BS) deployments. The viability 
of mmWave for energy harvesting though is unclear, due to 
the differences in propagation characteristics such as extreme 
sensitivity to building blockages. This paper considers a scenario 
where low-power devices extract energy and/or information from 
the mmWave signals. Using stochastic geometry, analytical ex¬ 
pressions are derived for the energy coverage probability, the av¬ 
erage harvested power, and the overall (energy-and-lnformation) 
coverage probability at a typical wireless-powered device in terms 
of the BS density, the antenna geometry parameters, and the 
channel parameters. Numerical results reveal several network 
and device level design insights. At the BSs, optimizing the 
antenna geometry parameters such as beamwidth can maximize 
the network-wide energy coverage for a given user population. At 
the device level, the performance can be substantially improved 
by optimally splitting the received signal for energy and infor¬ 
mation extraction, and by deploying multi-antenna arrays. For 
the latter, an efficient low-power multi-antenna mmWave receiver 
architecture is proposed for simultaneous energy and information 
transfer. Overall, simulation results suggest that mmWave energy 
harvesting generally outperforms lower frequency solutions. 

I. Introduction 

Millimeter wave (mmWave) communications is a key can¬ 
didate technology for future 5G cellular networks. This is 
mainly due to the availability of large spectrum resources at 
higher frequencies, which leads to much higher data rates. 
Recent research suggests that mmWave systems will typically 
feature (i) large-dimensional antenna arrays with directional 
beamforming at the transmitter/receiver—which is motivated 
by the small wavelength that allows packing a large number 
of antenna elements into small form-factors; and (ii) a dense 
deployment of base stations (BSs) to ensure a comparable 
coverage to ultra high frequency (UHF) networks & @- 
These mmWave design features are also attractive for RF 
(radio frequency) energy harvesting where a harvesting device 
may extract energy from the incident RF signals Q. This 
could potentially power the massive number of low-power 
wireless devices in future paradigms such as the Internet of 
Things 0. The signal propagation at mmWave frequencies, 
however, suffers from poor penetration and diffraction char¬ 
acteristics, making it sensitive to blockage by buildings |^, 
©■ It is, therefore, unclear if mmWave cellular networks will 
be more favorable for RF energy harvesting compared to the 
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conventional (below 6 GHz) frequencies. Further, the network 
level design principles for mmWave energy harvesting systems 
are not well understood. This motivates a network view of 
energy harvesting in a mmWave cellular network. 

A. Contributions 

In this paper, we provide a tractable framework to character¬ 
ize the performance of wireless energy and information trans¬ 
fer aided by a large-scale mmWave cellular network. Our anal¬ 
ysis accounts for the key distinguishing features of mmWave 
systems, namely the sensitivity to blockage and the use of 
potentially large antenna arrays at the transmitter/receiver. 
We first consider mmWave energy harvesting, where devices 
only extract energy from the incident mmWave signals. Our 
analysis models two operating scenarios, one where devices 
have their beams aligned to that of a mmWave BS, and 
the other where no such beam alignment is assumed. For 
both operating modes, we derive simple analytical expressions 
for metrics such as the energy coverage probability and the 
average harvested power using tools from stochastic geometry. 
We then extend the analysis to characterize the overall (energy- 
and-information) coverage probability for the general case 
where a device extracts both energy and information from the 
mmWave signals. 

To get design insights, we examine the network level perfor¬ 
mance trends in terms of key parameters such as the mmWave 
network density and the antenna geometry parameters for both 
operating modes of the energy harvesting devices. Numerical 
results suggest that narrower antenna beams are preferred 
when the users are aligned with a BS, whereas wider beams are 
favorable when no beam alignment is assumed. Our findings 
also suggest that there typically exists an optimum transmit 
antenna beamwidth that maximizes the network-wide energy 
coverage for a given user population. This implies that the 
mmWave BSs will need to adapt the antenna beam patterns 
depending on the fraction of the users operating in each mode. 

Similar to the BS-related parameters, we also investigate 
the role of the device-related parameters on the system per¬ 
formance. For example, the overall (energy-and-information) 
coverage probability can be improved by optimizing over a 
design parameter (power splitting ratio) to optimally portion 
the received signal between the energy harvesting and the in¬ 
formation decoding modules. Another important design feature 
at the user is the receive antenna array. Similar to the BSs, 
the mmWave users can, in principle, benefit from using large 
antenna arrays. For low-power energy harvesting devices, how¬ 
ever, the associated antenna circuity could increase the power 
consumption, offsetting the potential gains of large antenna 
arrays. To leverage multiple antennas at the receiver without 
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resorting to power-hungry circuit components, we propose 
a simple switch-based receiver architecture for simultaneous 
energy and information transfer. Simulation results reveal 
that the proposed low-power solution performs reasonably 
well compared to more advanced but power-hungry receiver 
architectures. 


The paper is organized as follows. In Section we in¬ 
troduce the system model. Section m presents the analytical 
results for mmWave energy harvesting. The case with simul¬ 
taneous information and energy transfer is treated in Section 
[Tv] We conclude the paper in Section |V] 

II. System Model 


B. Related Work 

Wireless energy harvesting is becoming increasingly fea¬ 
sible due to the reduction in the power consumption re¬ 
quirements of wireless sensors and the improvements in 
energy harvesting technologies 0-G3- This has also led 
to considerable research in advancing the theoretical under¬ 
standing of wireless-powered systems (see Q, pT] for a 
comprehensive overview). For example, wireless energy and 
information transfer has been studied for different information- 
theoretic setups such as a broadcast channel GD, a fading 
channel p3| , and an interference channel 0- Many of these 
papers highlight the fundamental trade-off between energy and 
information transfer efficiency and characterize the achievable 
rate-energy regions for different practical receiver architectures 

HD- 

Wireless energy and/or information transfer in large-scale 
networks has also been investigated |[T5)-pO). In p3) , the 
performance of ambient RF energy harvesting was charac¬ 
terized using tools from stochastic geometry. Using a repul¬ 
sive point process to model RF transmitters, it was shown 
that more repulsion helps improve the performance at an 
energy harvester for a given transmitter density. In | [T6) , p7) , 
cognitive radio networks were considered, and opportunistic 
wireless energy harvesting was proposed and analyzed. In 
a hybrid cellular network architecture was proposed to enable 
wireless power transfer for mobiles. In particular, an uplink 
cellular network was overlaid with power beacons and trade¬ 
offs between the transmit power and deployment densities 
were investigated under an outage constraint on the data links. 
A broadband wireless network with transmit beamforming was 


In this section, we introduce the network and channel 
models, followed by a description of the antenna model. The 
parameters defined in this section are summarized in Table 

A. Network Model 

We consider a large-scale cellular network consisting of 
mmWave BSs and a population of wireless-powered devices 
(or users) that operate by extracting energy and/or information 
in the mmWave band. The mmWave BSs are located according 
to a homogeneous Poisson point process (PPP) <I>(A) of density 
A. The user population is drawn from another homogeneous 
PPP $u(Au) of density A^, independently of $. In general, 
mmWave BSs and users may be located outdoors or indoors. 
Empirical evidence suggests that mmWave signals exhibit high 
penetration losses for many common building materials |^, 
|21|. Assuming the building blockages to be impenetrable, 
we focus on the case where the BSs and users are located 
outdoors. We say that a BS-user link is line-of-sight (LOS) 
or non-line-of-sight (NLOS) depending on whether or not 
it is intersected by a building blockage. Channel measure¬ 
ment campaigns have reported markedly different propaga¬ 
tion characteristics for LOS/NLOS links Q, To model 


blockage due to buildings, we leverage the results in 1231 


considered in |19|, where optimal power control algorithms 
were devised for improving the throughput and power transfer 
efficiency. Simultaneous information and energy transfer in a 
relay-aided network was considered in pO) . Under a random 
relay selection strategy, the network-level performance was 
characterized in terms of the relay density and the relay 
selection area. 

Our work differs from the prior work in that we inves¬ 
tigate wireless energy and information transfer in a large- 
scale mmWave cellular network. Due to the different physical 
characteristics and design features at mmWave, prior work on 
energy/information transfer in lower frequency networks does 
not directly apply to mmWave networks. In another line of 
work, the performance of mmWave cellular networks in terms 
of signal-to-interference-and-noise ratio (SINR) coverage and 
rate has also been analyzed using stochastic geometry GD. 
| |2^ . None of this work on mmWave networks, however, 
provides a performance characterization from the perspective 
of wireless energy and information transfer. 


where the buildings are drawn from a boolean stochastic point 
process. We define a line-of-sight (LOS) probability function 
p{r) = e~^'" for a link of length r, where ^ is a constant 
that depends on the geometry and density of the building 
blockage process: a BS-receiver link of length r is declared 
LOS with a probability p{r), independently of other links. 
While conducting stochastic geometry analysis, we will apply 
this result to split the BS PPP into two independent but non- 
homogeneous PPPs consisting of LOS and NLOS BSs. 

We allow the user population to consist of two types of 
users, namely connected and nonconnected. A connected user 
is assumed to be tagged with the BS, either LOS or NLOS, that 
maximizes the average received power at that user. Moreover, 
for the connected case, we assume perfect beam alignment 
between a BS and its tagged user, i.e., the BS and user 
point their beams so as to have the maximum directivity gain. 
Further, we assume that a BS serves only one connected user 
at a given time. For a nonconnected user, we do not assume 
any prior beam alignment with a BS, i.e., it is not tagged with 
any BS. This allows us to model a wide range of scenarios. 
For instance, due to limited resources, the mmWave network 
may (directly) serve only a fraction of the user population as 
connected users, leaving the rest in the nonconnected mode. 
Another interpretation could be that due to the challenges 
associated with channel acquisition, not all the users could 
be simultaneously served in the connected mode. We let e be 
the probability that a randomly selected node is a connected 
user, independently of other nodes. With this assumption, we 
can thin the user PPP into two independent PPPs $ii,con 
and $u,ncon, with respective densities eA^ and (1 — €)Xu- 





Note that an arbitrary user, either connected or nonconnected, 
may experience an energy outage if the received power falls 
short of a required threshold ip. This threshold would depend 
on the power consumption requirements of the receiver. To 
capture the sensitivity requirements of the harvesting circuit, 
we define ^min to be the harvester activation threshold, i.e., the 
minimum received energy needed to activate the harvesting 
circuit (the energy outage threshold ip would typically be 
greater than ipmin)- We use ^ to denote the rectifier efficiency. 
We define Peon (A, ipcoa) to be the energy coverage probability 
given an outage threshold '0con for a connected user, while 
Pneon (A, ipncon) denotes the same for the nonconnected case. 
With these definitions, we can define the overall energy 
coverage probability A(e, A,'0con, V^ncon) of the network as 

A{e,X,ip con 1 V^ncon ) = ePcon (A, ■i/’con) + (1 “ e)P neon {X,ip neon ) (1) 


where the energy coverage probability is a function of several 
parameters such as the BS density, the channel propagation 
parameters, as well as the antenna beam patterns at the trans¬ 
mitter/receiver. For cleaner exposition, we drop the subscript in 
'ipcon or Ipncon, using the notation A(e, A, ip) when the context 
is clear. In Section III-A we provide analytical expressions 
to compute the energy coverage probability in a mm Wave 
network. 


B. Channel Model 


We now describe the channel model for an arbitrary user 
without losing generality. Empirical evidence suggests that 
mmWave frequencies exhibit different propagation character¬ 
istics for the LOS/NLOS links ||^. While the LOS mm Wave 
signals propagate as if in free space, the NLOS mmWave 
signals typically exhibit a higher path loss exponent (and 
additional shadowing) 0 - We let au and on be the path 
loss exponents for the LOS and NLOS links respectively. We 
define the distance-dependent path loss for a user located a 
distance ri from the ^-th BS; gi{ri) = when the 

link is LOS, where the constant Cl is the path loss intercept; 
and gt{ri) = Cn^^^” for the NLOS case. Note that by 


including blockages in our model (Section II-A i, we capture 
the distance-dependent signal attenuation due to buildings. 
To simplify the analysis, we do not include additional forms 
of shadowing in our model. We further define hg to be the 
small-scale fading coefficient corresponding to a BS £ S $. 
Assuming independent Nakagami fading for each link, the 
small-scale fading power = \he\'^ can be modeled as a 
normalized Gamma random variable, i.e., ~r (AA, 1 /A^l) 
when the link is LOS and ~r (A^n, 1/A^n) for the NLOS 
case, where the fading parameters Wl and A^n are assumed 
to be integers for simplicity. 
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Fig. 1. Sectored antenna model. The antenna beam pattern is parameterized 
by the directivity gains for the main lobe (M) and side lobe (m), and the 
half power beamwidths for the main lobe (6) and side lobe (9). 


pattern, where (p gives the angle from the Foresight direction, 
M denotes the directivity gain and 9 the half power beamwidth 
for the main lobe, while m and 9 give the corresponding 
parameters for the side lobe. With this notation, mi e, A (’) 
denotes the antenna beam pattern at an arbitrary BS in <!>, and 
Am^ m, e, eA') denotes the same for an energy harvesting user 
in We further define Se = (<^f)AMi.mi.ei.e, (</'f )> 

the total directivity gain for the link between the £-th BS and 
the typical user; (pf and (p^ give the angle-of-arrival and angle- 
of-departure of the signal. 

Without any further assumptions about the beam align¬ 
ment between a user and its BS, we model the directiv¬ 
ity gain 5i, as a random variable. We assume the angles 
(pl and (p^ are uniformly distributed in [0,27r). Due to the 
sectored antenna model, the random variable 5^ = Di 
with a probability pi {i G {1,2,3,4,5}), where G 
{MtMr, Mtmi-, mtMr, mtmr, 0} with corresponding probabil¬ 
ities Pi G {qtqr,gtgT,qtqr,qtqr,qo}; the constants = ^, 
qt = ^,qT= = 2-gt-gt-<?r-gr- Note 

that £>5 = 0 models the extreme case where the BS and user 
beams have no alignment at all. Note that for the connected 
mode, since we assume perfect beam alignment between the 
typical user and its serving BS (hereby denoted by subscript 
0), the directivity gain Sg = MtMy due to the sectored antenna 
model. 


III. MmWave with Energy Harvesting 


In this section, we assume that each user is equipped with an 
energy harvesting circuit, and attempts to extract energy from 
the incident mmWave signals. No decoding of information 
is considered in this section. The case with simultaneous 


information and power transfer is treated in Section IV We 
first provide analytical expressions to evaluate the energy 
coverage probabilities for both connected and nonconnected 
users. We then validate the analytical model, and conclude 
the section by providing network level design insights. 


C. Antenna Model 

To compensate for higher propagation losses, mmWave BSs 
will use large directional antennas arrays. We assume that the 
BSs and users are equipped with A) and Ny antenna elements 
each. To simplify the analysis while capturing the key antenna 
characteristics, we use the sectored antenna model of Lig. [T] 
(except for Section IV), similar to the one considered in plF 
1^. We use Aj^ yyy g g{(p) to characterize the antenna beam 


A. Stochastic Geometry Analysis 

We first provide some lemmas before stating the main 
analytical results for this section. 

Lemma 1 (Modified from ( [2i] Theorem 8]): The probabil¬ 
ity density function (PDL) of the distance from an en¬ 
ergy harvesting user to its nearest LOS BS, given that the 
user observes at least one LOS BS, is given by tl (x) = 
2TrXBi^~^xp{x)e~'^'^^-l^o where x > 0 and Hl = 1 — 





TABLE I 

Model Parameters 


Notation 

Description 

N„ N, 

Antenna array size at the transmitter 
(t) and receiver (r) 

Mt, M- 

Main lobe directivity gain 

nil, 

Side lobe directivity gain 

6l, 9y 

Main lobe half power beamwidth 

9i, 9, 

Side lobe half power beamwidth 

$(A) 

BS PPP with density A 


User PPP with density A„ 

e 

Fraction of connected users 


Energy outage threshold 

fimia 

Harvester activation threshold 


Rectifier efficiency 

A(e, A,V') 

Energy coverage probability 

p{r) 

LOS probability function 

P 

Building blockage parameter 

(Tl, On 

LOS/NLOS path loss exponents 

Cl, Cn 

LOS/NLOS path loss intercepts 


LOS/NLOS fading parameters 

Pt 

Transmit power of BSs in $ 


^-2v\fo vp(v)dv jjjg probability that the receiver observes 
at least one LOS BS. Similarly, the distance distribution of 
the link between the user and its nearest NLOS BS, given 
that the user observes at least one NLOS BS, is given by 
tn (x) = 27rAi3N~^a:(l — where 

a; > 0 and = 1 — is the probability 

that the user observes at least one NLOS BS. 

Lemma 2 (Modified from j |27| Lemma 2]): Let pL and pN 
denote the probability that the energy harvesting user is 
connected to a LOS and a NLOS BS respectively, then pL is 
given by pL = Sl /q°° \i-p{‘v))vdv^_^ (a:)da;, where 

Pl{x) = a;^ and pN = 1 - pl- 

Lemma 3 (Modified from j |27| Lemma 3]): Given that the 
energy harvesting user is connected to a LOS mmWave BS, the 
PDF of the link distance is given by the expression tl (a;) = 
gLT-L(x) - 2 TrXff^'- \i-p{v)}vdv y^jjgrg a; > 0. Given that the 

e-L ’ 

user is connected to a NLOS mmWave BS, the PDF of the 
link distance is given by tn (a:) = * p(v)vdv 

for a; > 0 and pi^{x) = x'^. 

Leveraging Slivnyak’s theorem p5) , we conduct the anal¬ 
ysis at a typical energy harvesting user located at the origin 
without losing generality. We let Pt be the BS transmit power, 
and Y = X]rG$(A) Pt^iHtgi(ri) be the power received at the 
user. Recall that '0niin denotes the harvester activation threshold 
defined in Section |II-A| The energy harvested at a typical 
receiver (in unit time) can be expressed as 

7 = (2) 

where f G (0,1] is the rectifier efficiency. Note that we have 


neglected the noise term since it is extremely small relative 
to the aggregate received signal. The remaining parameters 
follow from Section [n| Recall that given a BS 7 S $(A), the 
corresponding fading parameters will be distinct depending 
on whether the link is LOS or NLOS, which in turn depends 
on the LOS probability function (Section II-Ai. Further note 
that for the connected case, it follows from Section |II-C| that 
(5o = AftMr for the link from the serving BS (denoted by 
subscript 0). 


Connected case: The following theorem provides an 
analytical expression for the energy coverage probability 
Peon (A, 4>) = Pr {7 > ^/>} at a connected user, where the 
random variable 7 is given in Q, and is the energy outage 
threshold. Note that Peon (A, V') can also be interpreted as the 
complementary cumulative distribution function (CCDF) of 
the harvested energy. 


Theorem 1: In a mmWave network with density A, the 
energy coverage probability Peon (A, fi) for the connected case 
given an energy outage threshold can be evaluated as 

Peon (A, fi') — Pcon.L ^A, Pl T Pcon.N ^A, Pn, (3) 

where = max Pl = 1 — Pn is given in Lemma 

2, while Pcon.L (•) and Pcon.N (•) are the conditional energy 
coverage probabilities given the serving BS is LOS or NLOS. 
These terms can be tightly approximated as 

Pcon,L (A, ^ 

oo 

Cfc (r) dr, (4) 

’’s 

where (x) = ^1 -|- ^ , the approximation 

constant a = N{Nl)~^ where N denotes the number of 
terms in the approximation, while Vg defines the minimum 
link distance and is included to avoid unbounded path loss at 
the receiver. Similarly, 

N 

Pcon,N(A,V')«5I(-l)" 

k^O 
oo 

J Ck (r) dr, (5) 

’’s 

where = + 

4 / 

Tfe,i (A,i/’,a;) = SttA^P, J 




X p{t)tdt, 


( 6 ) 

































4 ^ / 

Tfe,2 (A, V’,®) = 27rA^Pi J 


1 + 


aPtkDiCti 


-JVn 


X {1 - p{t))tdt, (7) 


and the distance distributions fL(-) and %(•) follow from 
Lemma 3. 

Proof: See Appendix A. ■ 

Recall that p(f) = is the LOS probability function 

defined in Section II-A| and captures the effect of building 
blockages. In (j^, the term (•) models the contribution from 
the LOS serving link, i (•) accounts for the other LOS 
links, and 2 (•) captures the effect of the NLOS links. Note 
that the ith term in (|^, Q corresponds to the contributions 
from the BS-user links having directivity gain Di. Similarly, 
Cf’ (•) in ([^ models the case where the serving BS is NLOS. 
Note that these terms further depend on the channel propa¬ 
gation conditions (aL, ctN, -^n, C'l, C'n), the network 
density A as well as the antenna geometry parameters (via 
Di, Pi), which are summarized in Table |I] Furthermore, the 
outage threshold ip will depend on the power requirements 
at a particular user, and would typically be greater than the 
sensitivity of the harvesting circuit (i.e., ip > ip^in such 
that Ip = j). Though the expressions in Theorem 1 can be 
evaluated using numerical tools, this could be tedious due to 
the presence of multiple integrals. To address this, we simplify 
the analysis by approximating the LOS probability function 
with a step function, and by further ignoring the small-scale 
fading. This result in a much simpler expression for the energy 
coverage probability. 
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Proposition 1: Let Rb = 

A7ri?|e-^’^^B, and W^k = for i e {1,2,3,4}, 

k G jo, 1, • • • ,iV}. The energy coverage probability can be 
further approximated as 


-Peon (A,'0) 





,W,kR. 


dt, 


( 8 ) 


where r(h;u,v) = a;^“^e““da; is the generalized incom¬ 

plete Gamma function. 

Proof: See Appendix B. ■ 

We note that Proposition 1 is relatively efficient to compute 
as it involves integration over a finite interval only, and 
because Gamma function can be readily evaluated using most 
numerical tools. We also observe that the coverage is mainly 
influenced by the LOS BSs. For example, a key term in Q 
is XttR^ which represents the average number of LOS BSs 
seen by the user. We now provide analytical expressions for 
the average harvested power at a connected user. 

Proposition 2: The average harvested power for the con¬ 
nected case Peon (A, ip) for an energy outage threshold ip G 

[lpn,in,Oo) is given by Foon(X,1p) = Poon (X, x) dx + 


V’Peon (A, "0) . 

Proof: The proof follows by noting that 7 has nonnegative 
support, and by treating Peon (■, •) as the CCDF of 7. ■ 

Here, Peon (•, •) follows from Theorem 1 or Proposition 1. 
Peon (A, Ip) can be interpreted as the useful average harvested 
power. This is because only those incident signals that meet the 
activation threshold can be harvested. To get further insights, 
we now analyze the limiting case ip ^ 0 of Proposition 2. This 
provides an upper bound on the average harvested power. 

Corollary 1: The average harvested power fop the limiting 
case lim Peon (A, V') = -Peon (A, 0) = c (pl-Pl + ^n-Pn), 

■0-J-O 

where 

00 

Pl = y (PtMtM.C'Lr-"^ -h (r) + (pL(r))) R (r) dr, 

(9) 


Pn = y (PMtM,GNr-“” -h (pN(r)) -f (r)) % (r) dr, 

^9 

( 10 ) 

4 pQO 

jx) = kCl DiPx / ( 11 ) 

4 'n (x) = kGn ^ DiPi ( ^ -/ p(i)dt j , 

( 12 ) 

and K = 27rAPt. 

Proof: See Appendix C. ■ 

and Pn denote the average harvested power given the 
user is tagged to an LOS or an NLOS BS. Note that the 
average harvested power is independent of the small-scale 
fading parameters. To reveal further insights, we provide 
the following approximation for the average harvested power 
(which is validated in Section [Ill-B| i. 

Corollary 2: The average harvested power for the limiting 
case. Peon (A,0), can be further approximated as 

Rb 2 

(a) f 

Peon (A, 0) « kM.M.Cl / 

’’s 

r (1 — O.Sol; Ayrrg, 00) — F (l — O.Sol; Xt:R%, 00) 

~ 2 {K.MtMrCL)~^ (A7r)i-0-5“L ■ 

Proof: The proof follows by using the simplifying as¬ 
sumptions of Appendix B, and by further ignoring the contri¬ 
butions from all but the serving LOS BS. ■ 

This approximation suggests that the average harvested power 
is mainly determined by the LOS serving link. Note that 
Pcon{‘,') grows linearly with the transmit power Pt since 
K = 27rAPt. Depending on the path loss exponent a^, it may 
exhibit a sublinear to approximately-linear scaling with the 
BS density A. When is large, the denominator 7“^" ^ in 
(a) overshadows the impact of A on the numerator . 

Therefore, the scaling behavior is essentially determined by 















K = 2tt\Pi, which is linear in A. This suggests that increasing 
the transmit power or BS density has almost the same effect 
on the average harvested power when ol is large (e.g., when 
ttL = 3). Also note that 0 is relatively simple as it is 
expressed in terms of the incomplete Gamma function only. 

Nonconnected case: Having discussed the connected case, 
we now consider the case where a user operates in the noncon¬ 
nected mode. The following theorem characterizes the energy 
coverage probability at a typical user for the nonconnected 
case. 

Theorem 2: In a mmWave network of density A, the energy 
coverage probability for the nonconnected case Pncon (A, ip) 
given an outage threshold ip can be evaluated using 

N 

P„co„(A,r/;)«^(-l)" 
k^O 

(14) 

where Tk,i (•) and 2 (•) are given by ^ and Q respec¬ 
tively, Ip = max , and Vg is the minimum link 

distance. 

Proof: The proof follows from that of Theorem 1 and is 
therefore omitted. ■ 

Similar to the connected case, the energy coverage probability 
for this case is also a function of the propagation conditions, 
the network density and the antenna geometry parameters. 
We note that the expressions in Theorem 2 are efficient to 
compute, obviating the need for further simplification. We now 
consider the average harvested power for the nonconnected 
case. 

Proposition 3: The average harvested power for the non¬ 
connected case Pncon (A, Ip) for an energy outage threshold 
V’ ^ [^min, 00 ) is given by P^con (A, Pncon (A, x) dx + 

V’-^ncon (A, . 

Proof: The proof follows from that of Proposition 2 and 
is therefore omitted. ■ 




Corollary 3: The average harvested power for the limiting 

case lim Pncon (A, Ip) = Pncon (A, 0) is given by 
-0—>-0 


(A,0) =C(«'l (rg)-f T-n (rg)) 


(15) 


where 'I'l (■) and T'n (•) are given in (111 and ( 12 i respec¬ 
tively. 

Proof: The proof follows from that of Corollary 1 and is 
therefore omitted. ■ 


The average harvested power for the nonconnected case scales 
linearly with the transmit power and the BS density. This 
follows from (15i as both 4 'l (•) and T'n (■) relate linearly 
with the transmit power and density via the term k = 27rAPt. 
This also suggests that increasing the transmit power or density 
has the same effect on the average harvested power. Note that 
this is different from the connected case where the path loss 
exponent affects how average harvested power scales with the 
BS density. 


B. Results and Design Insights 


We first verify the accuracy of the analytical expressions 


presented in Section III-A using simulations. We then study 
how key design parameters such as the antenna beam pattern 
affects the energy coverage probability in purely connected 
(e —> 1) and nonconnected (e —)• 0) networks. We also 
compare the performance of mmWave energy harvesting with 
lower frequency solutions. After developing key insights for 
purely connected/nonconnected scenarios, we provide energy 
coverage results for the general case (0 < e < 1), where the 
network serves both types of users. 

Validation: In the following plots, the users are assumed 
to be equipped with a single omnidirectional receive antenna, 
the mmWave carrier frequency is set to 28 GHz, the blockage 
constant /3 = 0.0071 ID, and Ip > ipinin- In other words, for 
a given ip, the plots are valid for any ipj^m < ip- Note that for 
the less relevant case when ip < ip^in, the energy coverage 
probability flattens out, and is specified by Peon (A, V'min) or 
Pncon (AjV'min)- Without loss of generality, we set the rectifier 
efficiency ^ = 1 since this parameter does not impact the 
shape of the results, i.e., setting ^ < 1 results in shifting all the 
curves to the left by the same amount. We assume the rectifier 
efficiency to be the same when comparing mmWave and UHF. 
Note that there are no standard values for ^ since prior work 
has reported widely varying values pO) , p6| depending on 
the device technology, operating frequency, etc. For example, 
ID suggest that a mmWave energy harvesting circuit 
may have better overall performance than its lower frequency 
counterparts. Fig. ^a\ and [2b| plot the energy coverage proba¬ 
bility for the connected case using different model parameters. 
The analytical results based on Theorem 1 are obtained using 
N = 5 terms in the approximation. The simulation results are 
generated using Monte Carlo simulations with 10,000 runs. 
Similarly, using Theorem 2, Fig. and plot the energy 
coverage probability for the nonconnected case. There is a 
nice agreement between analytical and simulation results. 


Connected case fe — 1): In Fig. 2a we plot the energy cov¬ 
erage probability with three distinct transmit beam patterns for 
a given network density. We observe that the energy harvest¬ 
ing performance improves with narrower beams, i.e., smaller 
beamwidths and larger directivity gains. As the beamwidth 
decreases, relatively fewer beams from the neighboring BSs 
would be incident on a typical user. But the beams that do 
reach, will have larger directivity gains, resulting in an overall 
performance improvement. This is possible due to the use of 
potentially large antenna arrays at the mmWave BSs. Note 
that this performance boost will possibly be limited due to 
the ensuing EIRP (equivalent isotropically radiated power) or 
other safety regulations on future mmWave systems | [28) . 

For the purpose of comparison, we also plot the energy 
coverage probability for UHF energy harvesting under realistic 
assumptions. Given the current state-of-the-art 0, the 
UHF BSs are assumed to have 8 transmit antennas each. 
Further, they are assumed to employ maximal ratio transmit 
beamforming to serve a connected user. For the channel model. 








we assume an IID Rayleigh fading environment and a path 
loss exponent of 3.6 (no blockage is considered). The network 
density is set to 25 nodes/km^, which corresponds to an 
average distance of about 113 m to the closest UHF BS. 
The carrier frequency is set to 2.1 GHz and the transmission 
bandwidth is 100 MHz. As can be seen from Fig. 2a mm Wave 
energy harvesting could provide considerable performance 
gain over its lower frequency counterpart. Moreover, the 
anticipated dense deployments of mmWave networks would 
further widen this gap. This effect is illustrated in Fig. 
where we plot the energy coverage probability for different 
mmWave network densities for a given transmit antenna beam 
pattern. In Fig. we use Proposition 2 to plot the average 
harvested power at a typical mmWave user against the transmit 
array size. The plots based on Corollary 1 are also included. 
We note that the limiting case ijj —i' 0 treated in Corollary 1 
closely approximates the average harvested power obtained 
using Proposition 2. This figure also confirms our earlier 
intuition that mmWave energy harvesting can benefit from 
(i) potentially large antenna arrays at the BSs, and (ii) high 
BS density, which would be the key ingredients of future 
mmWave cellular systems. Fig. 3b shows how the path loss 
exponent impacts the scaling behavior of the average harvested 
power with BS density, corroborating the discussion following 
Corollary 2. 

Nonconnected case (e —?► 0).' We now analyze the energy 
harvesting performance when the harvesting devices operate 
in the nonconnected mode. In a stark contrast to the connected 
case. Fig. shows that for the nonconnected case, mmWave 
energy harvesting could benefit from using wider beams. This 
is because BS connectivity (alignment) is critical for the 
nonconnected case. With wider beams, it is more likely that a 
mmWave BS gets aligned with a receiver, albeit at the expense 
of the beamforming gain. Furthermore, a comparison with 
UHF energy harvesting shows that mmWave energy harvesting 
gives a comparable performance to its UHF counterpart. 
Similarly, Fig. 4b plots the energy coverage probability for 
different deployment densities. We note that performance can 
be substantially improved with denser deployments, which 
would be a key feature of future mmWave cellular systems. 

General case (0 < e < 1).' Having presented the energy 
coverage trends for the two extreme network scenarios, we 
now consider the general case where the user population 
consists of both connected and nonconnected users. We expect 
this to be the likely scenario for reasons explained in the 
network model (Section [lI-Al l. As described in Section [ri-C| an 
antenna beam pattern can be characterized by the half power 
beamwidth and directivity gain for both the main and side 
lobes. By tuning these parameters, the beam pattern can be par¬ 
ticularized to a given antenna array. As an example, we assume 
that uniform linear arrays (ULA) are deployed at the mmWave 
BSs. We use the following relations to approximate the main 
and side lobe beamwidths as a function of the transmit array 
size: 6^ « — arcsin 




and 6lt « ^ 

^ TT 

We use Ml = lOUlog(Wt) and mi = V 


arcsin |^. 

(Ml — 12) foihme 


directivity gains of the main and side lobes pO) . To ensure 
the power normalization, the constant V is chosen to satisfy 
Am. + ^mi = 1. 

In Fig. 1^ we plot the overall energy coverage probability 
A(e,'(/),A) against transmit array size Ni for different values 
of parameter e. We find that the optimal transmit array size 
depends on the statistics of the user population. For example, 
when e is large, it is desirable to use large antenna arrays at 
the BSs. When e is small, it is favorable to use small antenna 
arrays to improve the overall energy coverage probability. 
Depending on the network load (or the user population mix) 
captured via e, the energy coverage probability can be sub¬ 
stantially improved by intelligent antenna switching schemes. 
Since the parameter e would typically vary over large time- 
scales, such schemes would be practically feasible. 

Having presented the energy coverage trends for mmWave 
energy harvesting, we now consider the scenario where the 
user attempts to extract both energy and information from the 
incident mmWave signals. 

IV. MmWave Simultaneous Information and Power 
Transfer 

In this section, we consider the case where the energy 
harvesting device also attempts to decode information from the 
received signals, in what is known as simultaneous wireless 
information and power transfer (SWIPT) 0^ HD- We now 
assume that the energy harvesting receiver is also equipped 
with an information decoding circuit. We focus on the case 
where a given user is already aligned with its serving BS, 
i.e., e = 1 for this section. Further, we consider a power 
splitting receiver architecture fTT ) where the received signal 
is split using factors ^/L' and yU — v, v G [0,1]. A fraction 
yr — 1/ of received signal is available for energy harvesting, 
while the remaining signal is used for information decoding. 
With this notation, the signal-to-interference-plus-noise ratio 
(SINR) at a typical receiver can be expressed as SINR = 
; where S = PiMiM^HQg()(ro) denotes the useful 
signal power and I = E£>o,ee<s>(\)\B(r,) PiSiHege(ri) gives 
the aggregate interference power from the neighboring BSs. 
is the thermal noise power before splitting, while cr^ captures 
possible signal degradation after power splitting. Similarly, 
7 = (1-+ I + l{s+i+cr^>.,p^i^} denotes the re¬ 

ceived signal power fed to the energy harvester. Note that 
a user will be in outage if the harvested energy and/or the 
SINR fall below their respective thresholds. We now define 
Hiuc(A, T, Tp, v) — Pr [SINR > T, 7 > r/;] to be the probability 
of successful reception given the SINR outage threshold T, 
the energy outage threshold ip, and the power splitting ratio 
v. Extending the results from the previous sections, we now 
provide an analytical expression to characterize the system 
performance with SWIPT. 

A. Stochastic Geometry Analysis 

Before stating the main result of this section, we first 
provide a lemma for the SINR coverage probability at a 
mmWave receiver HD- 


















Fig. 2. (a) Energy coverage probability A(e,‘i/;, A) for different transmit antenna beam patterns parameterized by [Mi,mi,9i,9t] in a purely connected 

network (e = 1, A = 100/km^). The performance improves with naiTower beams for this case. Pt = 13 dB, W = 100 MHz, ql = 2, = 4, A^l = 2, 

= 3, and Vg = 1 m. There is a nice agreement between Monte Carlo simulation (sim) results and the analytical (anlt) results obtained using Theorem 
1 with N = 6 terms, (b) Energy coverage probability A(l,‘i/;, A) for different network densities for connected users. Transmit beam pattern is fixed to 
[10,—10, 30°, 330°]. Other parameters are the same as given in Fig. Also included are the results based on the analytical approximation (approx) in 
Proposition 1. The approximation becomes tighter as the density increases. 


Lemma 4 (Modified from j |Z7] Theorem ]]): In a mm Wave 
network of density A, the SINR coverage probability 
Pcov (A, T, v) at a SWIPT device, given an SINR outage 
threshold T and a power splitting ratio v, is given by 

Pcov (A, r, v) = Pcov,L (A, T, v) Qi^ + Pcov,N (A, T, v) (16) 


Afc, 3 (r,a:) = 27rA^Pi f jl- 
Pn{ 2 :) 

X p{t)tdt, (21) 


1 + 


CNfcDiCLTa:“N 




where pl = 1 — Pn is defined in Lemma 2, and Pcov,l(') gives 
the conditional SINR coverage probability given the device is 
served by a LOS BS, and can be approximated as 




Pcov.L(A,r, 


k+1 
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oo 
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P Cl Mr 




(17) 


Similarly, the conditional SINR coverage probability for the 
NLOS case Pcov,n(-) is given by 


JVn 


Nfe+l 


Pcov,N(A,r,Z.)«y](-l) ^ ^ 

X (r) dr. 


AT \ r fcct^v°NT(<.^ + v 

' e PC-^MfMr 


(18) 


where 

Afc 


.1 {T,x) = 2t:X^p, j M 


1 + 


ci,kDiTx°‘^ 


iVLt“L 
X p(t)tdt, 


-JVl 


(19) 


Afe ,2 {T,x) = 2 -kX 
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ci,kDiCj^Tx°‘^ 
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X {l-p{t))tdt, (20) 


4 °° / 

Afc,4 (T,x) = 2-KX^pi J M 


1 + 


cyikDiTx°‘^ 


1 -iVN 


X (1 — p{t)) tdt, 


( 22 ) 


A = for i G {1,2, 3,4, 5}, cl = Wl (Wl!) and 

CN = IVn (7Vn!)“^. 

The following theorem provides the main analytical result of 
this section. 

Theorem 3: In a mmWave network of density A, the success 
probability P^a^{X,T,il},v) given the SINR outage threshold 
T, the energy outage threshold and the power splitting ratio 
v is given by 


Psuc (A, T, fi, v) Pcov (A, T, u) Peon (A, p) 

~t“ Peon (A.V?) [l Peon (A,m)], (23) 

where the SINR coverage probability Pcov(') can be evaluated 
using the expressions in Lemma while the energy coverage 
probability Pcon( ) follows from Theorem 1. We further dehne 

Pcon(A,/i) — Pcon,L (A,/i) Pl “P Pnn.NfA, where Tcon,L(’) 

and Pncon,N(’) are specified by (Wl and (m respectively, by 
setting Cfe{’) = C^(’) = f- Moreover, p and tp depend on 
several parameters including the power splitting ratio v, the 
SINR outage threshold T, the energy outage threshold ip, the 
harvester activation threshold ip^in, the rectiher efficiency f, 






























Fig. 3. (a) The average harvested power in a connected mmWave network for different number of BS antennas Nt and deployment densities A. Results 
based on Proposition 2 are obtained for -0 = —35 dBm. The analytical (anlt) results based on Corollary 1 ^ 0) are validated using Monte Carlo 

simulations (sim); and closely approximate the average harvested energy obtained using Proposition 2. The transmit antenna beam patterns ai‘e calculated 
using the approximations used for obtaining Fig.[^ Other simulation parameters are the same as used in Fig. For comparison, a plot for a UHF system 
is also included, (b) Plots the average harvested power (Corollary 1) vs. BS density for A^t = 32. The plot validates the approximation in jl3) . Moreover, it 
shows how the average power scales with the BS density for different path loss exponents ql. For illustration, also included are the solid lines corresponding 
to the (hypothetical) case when average power scales linearly with density. The scaling tends to become linear as ql is increased. 


and the noise parameters. Further, /i = — 

^ and i! = max 

Proof: See Appendix D. ■ 

Note that Pcon(A, p.) in (23 1 is the interference CCDF evaluated 
at parameter p. It plays a key role in determining the operating 
mode of the system. Though the interference is harmful 
for information decoding, it can be benehcial for energy 
harvesting. When the interference is high, the SINR coverage 
probability will typically limit the success probability. In the 
other extreme, the energy coverage probability will play the 
limiting role. Also note that the success probability can be 
optimized over the design paramter iz, given other parameters. 
Moreover, we can recover Theorem 1 and Lemma 4 from ( |2^ 
by letting ^ Q and T —)■ 0, respectively. 

Note that, in principle, the success probability at a connected 
mmWave energy harvesting or SWIFT device can be further 
improved by leveraging large antenna arrays at the receiver, 
thanks to smaller wavelengths. Though our analytical model 
allows the users to have receive antenna arrays, it implicitly 
assumes the presence of ideal RF combining circuitry con¬ 
sisting of power-hungry components such as phase shifters, 
multiple RF chains, etc. When large antenna arrays are used at 
the receiver, the power consumption due to additional antenna 
circuitry may get prohibitively high, overshadowing the array 
gains. As SWIFT typically targets low-power devices, we 
present a simple low-power receiver architecture in the next 
section. Note that the analytical results based on Theorem 3 
can be interpreted as an upper bound on performance when 
the receiver consists of suboptimal components (as is the case 
in the following section). 


B. Low-power Receiver Architecture 


We now propose a novel architecture for a mmWave SWIFT 
receiver with multi-antenna array, as depicted in Fig. § In 
this architecture, we assume per-antenna power splitting with 
parameter v (as dehned earlier). After power splitting, the 
input signal at each antenna passes through a rectiher, followed 
by a DC combiner that yields the harvested energy. For the 
information path, after passing through power splitters, the 
received signals are hrst combined in the RF domain using 
a combining vector w. The resulting signal is then decoded 
in the baseband. Because they require extremely small power, 
the combining vector is assumed to be implemented using 


switches j3l) , i.e., w = [wi, • • • , WjvJ* G [0,1]'^'. Note that 
both the signals for information decoding and energy harvest¬ 
ing are in the order of /rW (Fig. [7]). It is worth mentioning 
that recent results have shown that mmWave energy harvesting 
circuits can run with only a few pW | [2^ , |27|. 

We now derive the combining gain expression for the 
proposed SWIFT receiver architecture in Fig. Let y be 
the signal output at the RF combiner. If a BS applies a 
beamforming vector f G ([];tv,xi jq symbol s 

(where E [|sp] = Pfj to a target user, it follows that 


y = yfy [w*Hdfs -f w*rint + w*n] 


(24) 


where Hd G is the channel between the user and 

its serving BS, and Tint is the received signal due to the 
interfering BSs. Since the channel between each user and its 
BS is assumed to be single-path, the channel matrix Hd = 
ho^/go{ro)a.,. {fr) a* {ft), where at (ft) and {fr) are the 
array response vectors at the BS and user, respectively. Recall 
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Fig. 4. (a) Energy coverage probability A) for different transmit antenna beam patterns in a nonconnected network (e = 0, A = 100/km^). The 

performance improves with wider beams for this case. Other simulation parameters are same as given in Fig.[^ Monte Carlo simulation (sim) results validate 
the analytical (anlt) results obtained using Theorem 2 with N = 6 terms, (b) Energy coverage probability A(O,'0,A) for different network densities for 
nonconnected users. Transmit beam pattern is fixed to [10,—10, 30°, 330°]. Other parameters are same as given in Fig.[^ Monte Carlo simulation (sim) 
results validate the analytical (anlt) results obtained using Theorem 2 with N = b terms. 
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Fig. 5. The overall energy coverage probability A(e, "0, A) for different values 
of e. Depending on the fraction of users operating in connected/nonconnnected 
modes, the transmit aiTay size (which controls the beamforming beamwidth 
in this example) can be optimized to maximize the network-wide energy 
coverage. This could translate into massive gains given that the number 
of seiwed devices would be potentially large. The users are assumed to be 
equipped with a single omnidirectional receive antenna. The energy outage 
threshold -i/j is —70 dB for and —85 dB for T»u,ncon- P\ = 13 dB, 

A = 200/km^. Channel pai'ameters are the same as used in Fig. |2a| 



that go{ro) denotes the path gain from the serving BS, while 
(/)r and (pt respectively denote the channel angle-of-arrival and 
angle-of-departure at the user and BS. If the channel is known 
at the BS, and given the antenna model in Section |II-C| the 
BS will design the beamforming vector f to maximize the 
beamforming gain, i.e., to have |a* ((/)t)f|^ = Denoting 

a = /io\/(7o(^o)a-t ('/'t) f. the received signal in can be 



Fig. 6. Low power receiver architecture for SWIFT. 


written as 


y = y/v (aw*ar ((/>r) s + w*ri„t + w*n). 
The post-combining SINR can then be expressed as 


SINK = 


vPi\a\^ Nt |w*ai. {(j)r)\^ 

I + ’ 


(25) 


(26) 


where |w*ai. ((/)r)| represents the combining gain at the 
receiver, and I denotes the aggregate interference power. The 
SINR in p6| ) can be maximized if the receiver designs the 
optimum combining vector, which can be implemented by ac¬ 
tivating certain antennas on or off. This requires the receiver to 
have global channel knowledge, which is often challenging in 
practice. We relax this condition by assuming that the receiver 
has the angle-of-arrival information for the serving BS only. 
Ignoring the interference, we propose to design the combining 


vector by maximizing the SNR = 


Pi|a| A''t|w*ar(0.j,J| 


instead. 



































































Algorithm 1 Greedy Switch Combining Design 

Input Ar, 

Initialization w = 0, wi = 1 
for i = 2, • • • , Ar do 


if i 


^*- 1 ^ yy^gjfcd(n-l) cos( 0 ,) _|_ gjfcd(i-l) cos((^r) 


1 

i-1 


^^^]kd(n-l) cos{<l>r) 

Wj = 1 


2 

then 


2 


end if 
end for 


> 


i.e., the receiver designs its combining vector w such that 


w = arg max 
wG [0,1]^ 


|w*ai. ((/)r)r 


(27) 


The optimal solution to can be found by an exhaus¬ 
tive search over all possible combinations of w. For large 
receive antenna arrays, this could entail high computational 
costs, which would further increase the power consumption. 
Therefore, it is important to consider computationally efficient 
approaches for designing the combining vector. As outlined in 
Algorithm we propose a greedy solution for designing w 
by (step-wise) activating only those antennas that boost the re¬ 
ceived SNR. With w denoting the combining vector designed 
using Algorithm 1, the combining gain for the switch-based 

architecture can be dehned as Mr = -—— -m’)-^ 

|W|- 

where k denotes the wavenumber and a is the antenna element 
spacing. Despite its low-complexity, numerical simulations in 
the next section reveal that our low-power greedy approach 
could give a good combining gain, without losing substantial 
performance compared to more advanced but power-hungry 
solutions. 


C. Results 

Fig. 0 plots the overall success probability for a given 
transmit antenna beam pattern. The users are equipped with 
a single-antenna receiver, similar to the one in Fig. with 
Ar = 1. First, Fig. |7] shows that a reasonable success 
probability can be obtained with mmWave SWIFT system for 
typical mmWave propagation and system parameters. Further, 
this plot illustrates that the power splitting ratio v needs to be 
optimized for a given SINR outage threshold to maximize the 
overall success probability. Matching the intuition, the hgure 
shows that in the low SINR outage regime (when T is large), it 
is desirable to divert more power to the information decoding 
module, while a larger fraction of power needs to be portioned 
for the energy harvesting system in the high SINR outage 
regime (when T is small). This trend is consistent with prior 
studies on SWIFT architectures Ol- 

We now evaluate the performance of the proposed low- 
power receiver architecture for different number of receive 
antennas. In Fig. the success probability Psuc{^TT,tp,iy) 
is plotted for a hxed transmit antenna beam pattern. For the 
proposed architecture, the combining vector is obtained using 
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Fig. 7. A 3D plot showing the interplay between the success probability, 
the power splitting ratio u, and the SINR outage threshold T for a given 
energy outage threshold ip and network density A. As T gets large, the system 
becomes SINR-limited, and the optimum value of u increases, suggesting that 
a larger fraction of received signal should be used for information extraction to 
optimize the overall success probability. The transmit antenna beam pattern is 
set to yli 5 _i 5 iQcj 35 QO. Other parameters include i/i = — 70 dB, = —80 
dB, A = 200/km^, and Pt = 43 dBm. 


Algorithm 1, and the curves are averaged over the angle-of- 
arrival parameter. For comparison, we also plot the success 
probability for (fully digital) maximal ratio combining (MRC) 
receivers. We observe that the success probability improves 
with the receive antenna array size. Further, when the SINR 
outage threshold T is small, the success probability is mainly 
limited by the energy outage. This also explains why the 
success probability converges to a limit (determined by the 
energy outage threshold) as T decreases. Moreover, there are 
diminishing returns as the number of antennas are increased. 
A comparison with power-hungry MRC receivers shows that 
the proposed switch-based architecture performs reasonably 
well. This is particularly desirable for future mmWave SWIFT 
devices. 


V. Conclusions 

In this paper, we analyzed the energy harvesting perfor¬ 
mance at low-power devices powered by a mmWave cel¬ 
lular network. Using a stochastic geometry framework, we 
derived analytical expressions characterizing the performance 
of mmWave energy and information transfer in terms of 
system, channel and network parameters. Simulations results 
were used to validate the accuracy of the derived expressions. 
Leveraging the analytical framework, we also provided useful 
network and device level design insights. For the connected 
case when the transmitter and receiver beams are aligned, 
results show that the energy coverage improves with narrower 
beams. In contrast, wider beams provide better energy cov¬ 
erage when the receivers are not aligned with a particular 
transmitter. This trade-off is evident in the more general 
scenario having both types of receivers, where there typically 
exists an optimal beamforming beamwidth that maximizes 
the network-wide energy coverage. Moreover, we found that 













Fig. 8. The success probability for different number of receive antennas 
TVr at the user given a fixed transmit beam pattern ^ 15 , — 15 , 10 ° , 350 ° the 
BSs. Proposed low-power architecture achieves good performance compai'ed 
to superior receiver architectures. Other parameters include v = 0.5, ”0 = 
-70 dB, 0-2 = -80 dB, A = 200/km2, and P, = 43 dBm. 


several device-related parameters can significantly impact the 
system performance. For example, the performance can be 
substantially improved by optimizing over the power splitting 
ratio and by leveraging large antenna arrays. To allow using 
multiple antennas at the mmWave receivers while keeping the 
power consumption low, we proposed a low-power receiver 
architecture for mmWave energy and information transfer 
using antenna switches. Simulation results show that the 
proposed architecture can provide good gains for the overall 
mmWave energy harvesting performance. Simulation results 
also reveal that mmWave cellular networks could potentially 
provide better energy coverage than lower frequency solutions. 


Appendix A: Theorem 1 


The following inequality approximates the tail probability 
of a normalized Gamma distribution. 

Lemma 5 (From For a normalized Gamma random 

variable u with parameter N, the probability Pr {u < x) can 
be tightly upper-bounded by Pr (u < x) < (1 — 
where the constant a; > 0 and a = N{Nl)~«. 

We write Peon (A, i/j) = Pr 


Pr 


S + I>i; 


y > max = 

, where S = PiMiM^HQgQ(rQ) is the re¬ 
power from the serving BS, and I = 


ceived signal 

Sr>o,^G<E>(A)\B(rg) is the received signal power 

from all the other BSs. We can derive the result in Theorem 
by finding the conditional distributions Pcon.L (A, '0) and 
^con.N (A, V')- To proceed, first consider the conditional dis¬ 
tribution Pcon.L (A,?/') = Pr(S'-|- / >'0|L) given the receiver 
is aligned with a LOS BS (which is indicated by the subscript 
L in the following notation). 


Pcon,L (A, 'tp( — ^^5 j|L 


Pr u < 
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s,/lL (28) 


where we have included a dummy random variable 
in the first equation. Note that u converges 
to 1 as V —7> (X). Therefore, this substitution is in fact an 
approximation when N is finite. The introduction of u allows 
leveraging the inequality in Lemma 5, which leads to (a), 
where the constant a = N(Nl)~". The last equation follows 
from the Binomial series expansion of (6), and by further 
substituting a = To evaluate the expectation in (28i, 
consider 


Es,/|l 


g-Hs+i) 


Es|l [e-“''E,|s.L [e-“^]] • (29) 


The inner expectation in ( [29| can be simplified by applying 
the thinning theorem for a PPP | |25) . Note that $ can be 
independently thinned into two PPPs and <1)^, where the 
former comprises the LOS BSs whereas the latter consists 
of NLOS BSs. Therefore, we can interpret and as 
two independent tiers of BSs. The user will be tagged with 
either the closest BS in or in whichever maximizes the 
average received power at the user. We can further thin into 
four independent PPPs {‘hL}i=i, where each resulting PPP <1 >l 
contains BSs that correspond to a nonzero directivity gain Di 
with Pi being the thinning probability. This follows because the 
beam orientations are assumed to be independent across links. 
Thus, a link can have a directivity gain of Di with probability 
Pi independently of other links. We let the received power due 
to the transmission from the BSs in be Likewise, 
can be split into "'^h the corresponding received 

powers denoted by {7^}’^]^. Since the resulting PPPs are 
independent, (|29|) can be simplified as 


i=i 


E/|s,l [e-"'] = n®9S.L (30) 

where 
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— 27rXpi J I 1 


1 + aPtDiCi,re-°‘kNi, ^ 

\ Nt \ 

I p{t)tdt 


l + aPi£)^CLt~“L , 


(31) 


where (a) follows by conditioning on the length Tq of the 
serving LOS link, and by further noting that small-scale fading 
is independent across links. Here, B (to) denotes a circular 
disc of radius Vg centered at the typical user, (b) is obtained 
by using the moment generating function of a normalized 
Gamma random variable, while the last equation follows by 
invoking the probability generating functional p5) of the PPP 
Substituting (31 1 in the first (left) product term of (301 
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incomplete Gamma function. The result in Proposition 1 is 
obtained by assuming ps = 0 , and by further noting that the 
distance distribution simplifies to tl (a;) = 
the LOS ball approximation. 


27r^g Xttx J.Q 


. (32) 

By substituting ( [3^ in the second (right) product term of ([30| ) 
yields Q. Usingme expressions in (|30l l-((32]l in ( p9| l, ana by 
further evaluating the expectation of the resulting expression 
with respect to S, we obtain 


Appendix C: Corollary 1 

_ We derive Corollary_ 1 by finding the conditional means 
Pl = E [S' + /|L] and Pn = E [S + /|N]. Pl can be evaluated 
by conditioning on the link distance tq from the serving BS 
as follows. 


1 + aP,M,MrCLr-“L)VL“\ 

where we have again used definition of the moment generating 
function of a normalized Gamma distribution. i (•) and 
Tfc, 2 ( ) are given in ^ and Q respectively, denotes 
the minimum link distance, while the distance distribution 
is provided in Lemma 3. Using ( [28] l and ( [29| , we can thus 
retrieve the expression in 0 . We can similarly derive the 
conditional distribution Pcon.N (A, "0) = Pr (S + / > V’IN) in 
Q for the NLOS case. 

Appendix B; Proposition 1 
To simplify the analysis, we approximate the LOS probabil¬ 
ity function p(r) by a step function p{r) = l{o<r<RB}> 
the BSs within a LOS ball of radius are marked LC)S with 
probability 1, while the rest as NLOS pTl . The radius Ps is 
chosen such that the LOS association probability pl remains 
the same (as in our original model). Using a step function 
for p{r), it follows from Lemma 2 that pl = 1 — 

and Rb = . Moreover, we neglect small scale 

fading for all the links except for the serving BS. We also 
ignore the NLOS signals in the analysis. This effectively 
leads to a scenario where the user receives signals from the 
BSs within the LOS ball only. Intuitively, this would be 
the likely scenario in sufficiently dense networks. We only 
list the key steps since the rest of the proof follows from 
Appendix A. Ignoring the NLOS signals, we approximate ( |30| ) 
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(35) 


where (o) is obtained by averaging over the fading distribution, 
followed by invoking Campbell’s theorem p5) , while (|35l) 
follows from the definitions of T'l and T'n provided in 


TB 


and (12 1 respectively. Taking expectation of E[S' + /|ro,L] 
with respect to Tg using Lemma 3 yields 0 . The expression 
for Pn is (lOl can be derived using similar steps. 


Appendix D; Theorem 3 

From 0 ^ it follows that the harvested energy 7 = 

Let Y = {1 - v){S + I + 0 -^), where S = 
PtMiM,Hogo{ro) and I = J2e>o,ee<s>{\)\M{r,) PtStHege{n) 
respectively denote the contributions from the serving 
and the interfering BSs. To find Pjuc (A, P, i/>, v) = 
Pr [SINK > T, 7 > i/)], consider 
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(34) 


Here, (a) follows by ignoring the small scale fading and invok¬ 
ing the probability generating functional p?) of PPP, (b) by a 
change of variables, and (c) by the definition of the generalized 


Pr [7 < p] 

~ Pcov(A, P, /A)P;on(A, p') -p Pcon(A, p') |^1 -Peon (A.p)] (36) 

where the expectation in (a) is with respect to the interference 
7, (b) is obtained by further conditioning on I to be greater (or 
smaller) than ajiarameter p which follows from the inequality 
T (^I + — I — a^. The approximation (or 

effectively an upperbound) in (c) results from dropping the 
conditions I > p or I < p while calculating the expectation. 




















































Finally, the SINR coverage probability P^a^{\^T,v) follows 

from Lemmaj^ and the energy coverage probability Peon (A, (/?) 

from Theorem Pcon(A,/r) is the interference CCDF evalu¬ 
ated at p. 
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